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Abstract 

The mathematical formulation of the model for molecular move- 
ment of single motor proteins driven by cyclic biochemical reactions 
in an aqueous environment leads to a drifted Brownian motion char- 
acterized by coupled diffusion equations. In this article, we introduce 
the basic notion for the continuous model and review some asymptotic 
solutions for the problem. (For the lattice model see 0, Stochas- 
tic, nonequilibrium thermodynamic interpretations of the mathemati- 
cal equations and their solutions are presented. Some relevant mathe- 
matics, mainly in the field of stochastic processes, are discussed. 

Keywords: augmented Huxley equation; probability circulation; en- 
tropy production; macromolecular mechanics; Markov process; nano- 
biochemistry; nonequilibrium steady-state; protein; singular perturba- 
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1 Introduction 

One of the fascinating aspects of protein molecules in the biological world 
is their ability to perform various, almost "magic-like" tasks |T^. A partic- 
ular class of proteins known as molecular motors can move linearly along 
its designated track against an external force by utilizing the biochemical 
energy source, adenosine triphosphate (ATP). In this manner, the motor 
proteins act as miniature engines converting chemical energy to mechani- 
cal work. Movement of single protein molecules inside a cell, however, has 



to experience thermal agitation from the aqueous environment in the cy- 
tosol. The movement is therefore a Brownian motion with drift (convective 
diffusion) |^]. 

Such movement provides the molecular basis for muscle contraction and 
various cellular transport processes |24, 25]. Motor protein kinesin is known 



to carry out intra-cellular vesicle transport along microtubules. Various 
polymerases are moving along their corresponding templates. All these pro- 
cesses are essential to a living cell. In a muscle cell, the motor protein is 
called myosin, and its designated track is called an actin filament. The 
actin filament has a periodic structure of ~36nm. Therefore, without loss 
of generality, we assume that a myosin molecule moves in a force field with 
periodic potential energy function U{x): U{x + L) = U{x) where L is 36nm 
for actin. 

Treating the center of mass of the motor protein as a Brownian mo- 
tion with the presence of a periodic energy potential, its movement can be 
modeled by the Smoluchowski equation [^] 



dt dx ' dx'^ dx \ j3 

where D and (3 are, respectively, the diffusion and frictional coefficients. 
F[x) = —dU{x)/dx is the force of the potential C/, representing the molec- 
ular interaction between the motor protein and its track. P{x, t) is the 
probability density function of the motor protein at position x for time t. 
The first equality in Eqn. (|^) is a continuity equation in which J is a prob- 
ability fiux. The first term on the right-hand-side is associated with the 
diffusion fiux according to Fick's law. The second term is due to the con- 
vection associated with an overdamped Newtonian motion: f3x = F{x). In 
fact, Eqn. (|l|) is mathematically equivalent to an overdamped Newtonian 
motion with a white random force f{t) representing the incessant collisions 
between the motor protein and the water molecules: /3x = F{x) + f{t) ||6^ . 

The probability density P{x,t), as the solution to Eqn. (|l]), gives the 
mean position of the motor 

{x){t) = / xP{x,t)dx. (2) 



2 



Moreover, its velocity is related to the flux J{x,t): 



d , ^ dP(x,t) , 

-{xm = / x^^dx 



dt J-oo dt 

<-°° dJ(x,t) , 

X — dx 

ox 



oo 



oo 



J{x,t)dx. (3) 

When the motion of a motor protein becomes steady and if we are only 
interested in the mean velocity, we need only consider the steady-state so- 
lution for X € [0, L] with periodic boundary conditions. With this setting, 
the steady-state velocity of the motor protein movement is ^; = LJ. There- 
fore studies on the steady-state movement focus on the flux J. To flx our 
terminology, we will refer to a stationary solution as a steady-state, but a 
stationary solution with zero flux as an equilibrium. 

A little mathematical analysis immediately shows that with U{0) = 
U{L), i.e., Jq F{x)dx = 0, the stationary solution of (0) allows only zero 
flux (J = 0). Therefore, in a one- dimensional periodic structure, there is no 
driving force to bias an inert Brownian particle to move in either direction. 
The simplest model given by Eqn. (||) fails to capture the essence of motor 
protein movement. 

The driving force for a motor protein comes from a very important bio- 
chemical reaction, occurring inside the protein, called ATP hydrolysis: 

ATP + H2O i± ADP + Pi 

9 

where H2O is water, ADP is adenosine diphosphate, and Pi is phosphate. 
Chemical reactions like this are well characterized by a two-state Markov 
process (or more generally, m discrete states) 

dPATP rjy ,0 



dt 
dPADP 
dt 



fPATP - gPADP (4) 



where non-negative / and g are rate constants for the reaction |35]. The 
deeply insightful work of Sir Andrew Huxley in 1957 was to introduce in- 
ternal conformational states to the Brownian particle and to couple the 
(bio) chemical reaction in (^) with the motor protein movement in (|^) ||2~ 



This leads to the following equation [H, 43 , nn, M , known as a coupled dif 



fusion system in mathematics , for the movement of a Brownian particle 
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with internal structures and dynamics: 



dPjx, 
df 
dPjx, 

df 



D 



d^P{x, 



dx^ 



f3- dx \ dx 



P{x, -) - fix)P{x, -) + g{x)P{x, +) 



dx'^ (3+ dx \ dx 



P{x, +) + f{x)P{x, -) - g{x)P{x, +), 



(5) 

where f{x) and g{x) are non- negative periodic functions. In terms of this 
augmented Huxley model, the motor protein can be either attached to (+) or 
detached from (— ) the filament with respective interaction energy functions 
f7+(x) and U^{x) (In the original Huxle model, U-{x) = 0). D± {P±) are 
the diffusion (frictional) coefficients of the motor protein in the attached and 
detached states, respectively. The attach-detach transition is coupled to the 
ATP hydrolysis. Therefore, in (|5|) the biochemical reaction is coupled to the 
movement of the motor protein. More importantly, note that in either the 
attached or detached state, there is no bias for the motor protein movement! 
However, the attach-detach transition driven by the ATP hydrolysis leads to 
a biased motion of the motor protein (J^O). The chemical energy in ATP 
is converted to the mechanical motion of the motor protein. This we shall 
show. 

At this point, it is fascinating to read the now classic work of Huxley 
on the theory of muscle contraction, which was written decades before the 
discovery of the motor protein molecule in its individual form. The Huxley 
model works as follows (p. 281 of ||2^). Initially, the myosin (M) and actin 
filament (A) are detached; M oscillates (fluctuates) back and forth about 
its equilibrium position O as a result of thermal agitation (with diffusion 
coefficient D_ for the Brownian motion). If A happens to be within the 
range of positions where /, the rate of association, is not zero, there is a 
chance that combination will take place (this event happens with a proba- 
bility characterized by the Markov rate process); when this has happened 
the tension in the elastic element {F, i.e., the molecular interaction between 
M and A) will be exerted on the actin thread by M (and, conversely, the 
actin will exert the same force with opposite sign on M). As one can see, 
Eqn. (H) is the mathematical formulation of this model described in words. 
Readers who check Huxley's paper will find several differences between his 
original equation and (|5[). These differences arise because we have formu- 
lated a microscopic model for single motor proteins while a model for muscle 
contraction has to deal with a large number of myosin molecules. It can be 
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shown that when stringing many motor proteins into a rigid chain, Huxley's 
original equation can be derived from Eqn. (P) ||2^ . 

Since the work of Huxley, there have been many investigations following 
and expanding the basic notion of the Huxley equation. Most notably are 
the work of Hill |22, p^] who provided the Huxley equation with a sound 



thermodynamic basis, and the work of Astumian and Bier, and Peskin et 
al. who arrived at (^) from the Langevin dynamics (stochastic differential 
equation) point of view IQ, There is also a body of literature on Brownian 
ratchet, whose basic movement is characterized by equations identical to (^) 
[^, 0, 42, 1^]]^ In a recent paper, Bentil ||5| applied the Huxley model in 
conjunction with Langevin dynamics to simulate single myosin experiments. 
Qian established a relationship between coupled diffusion like (S) 

and the circulation of a Markov process and its entropy production |2^, ^ . 

Eqn. is certainly an oversimplified model for any realistic biolog- 
ical system. However, it captures the essence of a theory which unifies 
microscopic motor protein movement and macroscopic muscle contraction; 
thereby it provides a concrete model for chemomechanical energy transduc- 
tion in living organism. Hence it deserves further detailed investigation as 
a topic in biophysics and physical chemistry. The mathematical treatment 
of (I) has been mainly in terms of differential equations. However, the na- 
ture of Brownian motion of the motor protein also calls for a treatment of 
@ in terms of stochastic processes. We will review some of the pertinent 
mathematics in Section 4. 



2 Mathematical Analyses of Several Limiting Cases 
of Augmented Huxley Equation 

While the augmented Huxley Eqn. (||) is difficult to solve in general 
due to the non-local (non-equilibrium) nature of the steady-state, particular 
limiting cases can be analyzed to gain insights into the theoretical model. 
In this section, we present some known and also some new results. 



^It is important here to qualify the term "Brownian ratchet" which does not involve 
any temperature gradience [M . It is an isothermal device in which useful work is derived 
from non-equilibrium fluctuations. It violates the detailed balance due to active pumping 



5 



2.1 Limit of rapid biochemical cycling 

One particularly interesting limiting case is when the biochemical reactions 
are rapid with respect to the diffusion. Analysis of this limiting case clearly 
demonstrates how the internal biochemical reaction can give rise to a uni- 
directional motor protein movement. Hence it demonstrates the validity 
of mathematical models for motor proteins in terms of coupled diffusion 
equations. 

Consider Eqn. (|5|), rapid biochemical reaction means we have conditional 
probabilities :0 

Pi-\x) = ^r^^, ^^^^ 



f{x)+g{x)' f{x)+g{x) 

and thus P{x) = -P(x, — ) + -P(x, +) satisfies 

dP(x\ 9^ — d — 

m = d^2iDix)Pix,t)) - ^(F(x)P(x)) (6) 

in which 

^ f{x)D+ + g{x)D. ^ f{x)F^{x)+g{x)F^{x) 

^^''^ f{x)+g{x) ' "^^"^^ f{x)+g{x) 

Note that Eqn (^) is similar to Eqn. (||) but with one crucial difference: The 
mean force function now satisfies Jq F{x)dx ^ even though both F±{x) 
satisfy Jq F±{x)dx = 0. The potential of the average of forces of periodic 
L potentials is in general not periodic. This indicates that the biochemical 
reaction (^p provides a drift for the motor protein movement in a periodic 
system. Finally, the transport flux can be obtained by solving @ 

~ lo e-"(^)dx Jq e<^)dx - (1 - e«(-^)) /(f e-<^)dx Jq e«(^')dx' ' 

where u{x) = - f^T{x)dx. Therefore, if f^T{x)dx = -u{L) > 0, $ > 0. 

^Readers who are familiar with the method of singular perturbation will identify this 
problem. Since we are seeking a nontrivial solution for a homogeneous equation, the 
solution is not unique. The functions F± has at least two zeros at which boundary layers 
might be expected. For more detail on this type of equations, see |39f|. 



6 



2.2 Limit of rapid diffusion 



Another hmiting case which has been nicely analyzed by Peskin et al [43| is 
when the diffusion is very rapid in comparison to the Markovian transition 
(Brownian ratchet): there are rapid equilibria for P{x, — ) and P{x, +). 

The mathematical problem is framed as follows. Let's consider the sta- 
tionary coupled diffusion: 

d^Pix,—) d / dU-(x) ^. A /„/ XT-./ \ / \ T-./ 

+ -r\ , -) - <f{x)P{x, -) - g{x)P{x, +)) = 



dx"^ dx V dx 

^^^%^ + ^[^^P{x,+))+e{f{x)Pix,-)-g{x)^^^^ = 0, 

(8) 

where /(x), g{x), U±{x) are periodic functions. When the regular perturba- 
tion parameter e = 0, this system of uncoupled diffusion has zero transport 
flux For small e, we therefore have ^ = e(p + ... and the asymptotics can 
be obtained by the method of regular perturbation Q. At this point it is 
important to notice the other flux, the circular flux 11 which moves the mo- 
tor protein forward in the (-l-)-state but moves it backward in the (— )-state. 
Therefore 11 does not contribute to the net transport but only generates 



heat. This type of flux is known as futile cycle in muscle biochemistry |45]. 
Perturbation calculations show that p^ ] 

lo 9ix)l^+ix)dx Jq f{x)n-ix)dx ^ ^^^^^ 
lo 9{x)fJ'+{x)dx + Jq f{x)n-{x)dx 

and 



$ = e dx[v+{x) - u_{x)] r dx'[f{x')^i-{x')-g{x')ii+{x')]+0{e^) 
Jo Jo 



(10) 



where 

p-U-{x) p-U+ix) 



^ Joe-U-(-)dx' j^e-u+i-)dx' 

^U-(x) ^U+(x) 



It is interesting to note that if f{x)fi-{x) — g{x)fi+{x) = V x, then the 
system is reversible and the steady-state is in fact an (thermal) equilibrium. 
In mathematical term, the Markov process is symmetric ^2|. In applied 
mathematics, the symmetricity leads to the Grasman-Matkowsky variational 



method |3C]. More recent progress can be found in 
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2.3 Limit of the original Huxley model 

In the original Huxley model, the interaction between the track and the 
motor in the detached state is assumed to be zero. Hence U- (x) = in Eqn. 
@. Furthermore, it is also generally accepted that Z^^ << i.e., the 
Brownian motion of the motor protein in the attached state is negligible. 
The solution of (^) when — > is a problem of singular perturbations 
Note that because of the periodic U{x), F{x) has zeros on [0,L]. 
Thus the singular perturbation problem has at least two linear turning points 



[p8| , |3q| . The reduced equation when = is: 

'^^^^-f{x)P{x,-)+g{x)P{x,+) = 

-^{F{x)P{x,+)) + f{x)P{x,-)-g{x)P{x,+) = 0, (11) 

in which we have set -D_ = /?+ = 1 for simplicity. We are particularly 
interested in finding a condition for the existence of a solution corresponding 
to unidirectional motion. In the steady-state, the total transport flux of the 
system is a constant :0 

F{x)P{x,+)-^^^^^<^ (12) 
ax 



Using Eqn. (O) and eliminating P{x,+) from Eq. (11) we then have 



fip(r —\ dPlr —1 

F{x) ;j ^ + g{x) \' ^ - F{x)f{x)P{x, -) = -%(x), 
ax^ ax 



(13) 



where the inhomogeneous term <^ on the rhs is to be determined by the nor- 
malization condition J^[P{x,-'r)+P{x, —)]dx = 1. The boundary conditions 
for Eq. (|l3|) again are periodic. 

Eqn. ( [l3D has singular points at the zeros of F{x). A simple local 
analysis shows that for nonzero $ and a physically meaningful P{x, — ) > 0, 
the solution to Eqn. (13) has to be nonanalytic at these singularities. This 



nonanalytic behavior, however, is expected to be obviated in an asymptotic 
study of the full equation (^) with small . 

^We use $ to denote the transport flux J+ + J_. As will be shown below there is 
another type of flux, circular and non-transport 11 = J+ — J_ , in these systems . 
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3 Entropy Production in Nonequilibrium Steady- 
State 



We now give a brief discussion of the nonequihbrium thermodynamics 
in terms of Eqn. (^). Hih ||2^, ^] has given an extensive account of this 
subject. We only discuss some recent developments in connection with the 
notion of entropy production |4C]. The concept of entropy production rate 



(e.p.r.) can be easily introduced, mathematically, in terms of Eqn. (Q). 
The validity of this novel thermodynamics of nonequilibrium steady-state 
(NESS), however, remains to be experimentally tested. For more discussion 



see 40, ^7^, B4, 56 



Associated with (||) is a functional A[P(x)] called the Helmholtz free 
energy in thermal physics []58| , which in units /c^T (T is temperature and 
/cfi is the Boltzmann constant) is defined as: 



A[P{x,t)]= U{x)P{x,t)dx+ P{x,t)logP{x,t)dx. (14) 
Jo Jo 

When P{x,t) changes with time according to Eqn. (^, the production rate 
of total entropy is the rate of decrease in A of the system,^ which can be 
computed: 

dA fL dPix) fLdPix) 

e.p.r. = — = — U(x) — ; — dx — / — ; — iosP(x)dx 

^ dt Jo dt 70 

= / {U{x) +logP{x))^J{x)dx 
Jo dx 

= £ (^{F{x) - log P{x)J J{x)dx 
1 

= J^{x)P'\x)dx >0. (15) 

D Jo 

The last step used the definition for J given in Eq. (|l]). This is the second- 
law of thermodynamics in terms of Smoluchowski equation, which is the 

*From thermodynamics stand point, a macromolecule is an isothermal system in con- 
tact with a thermal environment (i.e., aqueous solution) with temperature T. There is 
energy (chemical and heat), but no material, exchange between the system and its environ- 
ment (clamped ATP and ADP concentration in the aqueous solution, and heat bath). The 
system and its environment as a whole is an isolated system with a constant total energy; 
and this is approximately hold also for a sufficiently large heat bath. In mathematical 
terms, we have dA/dt = dE /dt — TdS/dt where E is the internal energy of the system and 
S is the entropy of the system. With respect to the system and the environment together 
as a whole dA = {dEtot ~ dEenv) - T{dStot - dSenv) = -TdStot + {TdSenv - dEe 
—TdStot- For an isolated system (microcanonical ensemble), dE^nv/dSenv ~ T 



-^env } 
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nonequilibrium counterpart of a canonical ensemble in statistical mechanics. 



Eqn. (15) can be generalized to calculating the entropy production rate 
(e.p.r.), as well as the heat dissipation rate, in a NESS in which Stot continues 
to increase [^H]. To see this, we note the entropy of the system is defined as: 



S[P{x,t)] = -[ P{x,t)logP{x,t)dx. (16) 
Jo 

Therefore in NESS, 

S = - FJdx + e.p.r = 0. (17) 
Jo 

where F = —dU/dx, and the first term on the right-hand-side is the heat 
dissipation rate (h.d.r.). Therefore, in a NESS, the e.p.r. is equal to the 
h.d.r. 

Recent work in mathematical physics on entropy production in nonequi- 
librium systems |l9| also focus on appropriately setting up the non-equilibrium 
steady-state with an external force and a thermostat simultaneously acting 
on a Hamiltonian system. The force supplies energy while the thermostat 
removes heat in order to keep the system in a steady-state with bounded en- 
ergy. This leads to a random dynamical system in which entropy production 



is cogently defined |60|. The Smoluchowski approach we adopt has a quite 
similar setting: a driving force due to chemical reaction (rather than me- 
chanical force) and an implicit thermostat: The Smoluchowski equation is a 
consequence of an overdamped Newtonian system with Maxwellian distribu- 
tion for the velocity of the particles [^]. In fact, the diffusion coefficient D 
and frictional coefficient f3 in (|^) define the temperature of the thermostat: 
T = PD/ks- The interesting mathematical questions are when these ran- 
dom dynamical systems become diffusion processes and whether the entropy 
production proposed in these studies is equivalent to Eqn. ( [l5| ) for diffusion 
processes [^]. The recent work by Lebowitz and Spohn has provided 
some insights on this problem [^2|. In a different approach to weak random 
perturbation of Hamiltonian systems in a plane, Freidlin and Wentzell map 
the system to a diffusion process on a graph, which consists of vertices corre- 
sponding to the stationary states and edges corresponding to energy basins 
[jTsf. However, this approach remains to be generalized to higher dimensional 



Hamiltonian systems, and its relationship to the (Kramers') transition-state 



rate theory in theoretical chemistry |21| also remains to be elucidated 



It is interesting to note that in a steady-state, the flux J is a constant 
and all the entropy produced will become the dissipated heat. Hence e.p.r. = 
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(JV-D) Jq P''^{x)dx > J^L^/D which equals /JJ^L^ due to the Einsteine 
relation in ksT units D(3 = 1. The j3J'^L'^ term is the energy dissipation 
due to a deterministic motion with velocity v = LJ and frictional coefficient 
P in continuous medium. The inequality indicates the additional dissipation 
due to random motion. It also indicates that when P{x) = constant, the 
e.p.r is at its minimum, i.e., the chemomechanical energy transduction is at 
its maximal efficiency. 

To generalize the concept of entropy production to Eqn. (^) is math- 
ematically straightforward. This yields a novel thermodynamic theory for 
NESS, which is particularly relevant to motor proteins. The importance of 
the theory is that it relates e.p.r. to the heat production of a working motor, 
which is a quantity that can be experimentally measured. With some simple 
algebra, we have 



h.d.r. 



^ [f+(x)J+(x) + F_(x)J_(x) +i(x)log f44 

L \g{x) 



dx 

(18) 



rL r 



e.p.r. 



in which 



j2(x)P(x,+) + J!(x)P(x,-) + j(x)log 
dP{x,±) 



f{x)P{x,-y 
gix)P{x,+), 



dx 
(19) 



J±(x) 



dx 



+ F±{x)P{x,±), 



and 



j{x) = f{x)P{x,-) - g{x)P{x,+). 

It is obvious that e.p.r. is non-negative and equal to zero if and only if when 
detailed balance is hold |26, 4£]. In NESS without external load, h.d.r. = 
e.p.r. 

If there is an external load F^xt, then Eqn. ( p!8| ) can be further broken 
down into 



\[ADP][Pi] 



j{x)dx-Fext{J+ + J-)L. (20) 

JO 



Hence, The free energy from ATP hydrolysis is equal to exactly the sum of 
the work done against the external load, F^xtV-, and positive heat dissipation, 
Eqn. dH). 
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4 Some Relevant Mathematics on Coupled Diffu- 
sion 



While P{x,t) in Eqn. (||) characterizes a stochastic process in terms 
of probabihty density at each time t, P{x,t)dx = Prob{x < Xf < x + dx}, 
there is an alternative view of a stochastic process in terms of its trajectories. 
In this approach, all possible trajectories {Xf|t > 0} form a function space 
Q and a probability density (a measure) is defined. This naturally leads to 
the notion of "propagator" (a semi-group) which is formally defined as: 

P(x,t) = P(x,0)e^* (21) 

where the exponential operator acts on the distribution P{x, 0) as a "row 
vector" . The operator £ satisfies the backward Kolmogorov equation: 

^ = P(..t)^. (22) 

or its conjugate £* satisfies the forward Kolmogorov (or Fokker-Planck) 
equation: 



dP{x,t) 



£*P(x,t). (23) 



For symmetric operator £* = £. The symbolic relationship between the 
operator £ and the propagator e"^* has been made rigorous in terms of lin- 
ear operators in a Banach space and is known now as Hille-Yosida theorem. 
Hence the modern theory of Brownian motion has brought several math- 



ematical disciplines to bear |15, 29|: partial differential equations, linear 



operators on functional spaces, and harmonic analysis. 
4.1 Feynman-Kac formula 

A major result in this area is a relationship between the solution of a bound- 
ary value problem (BVP) and the mean exit time (first passage time) of a 
diffusion process. This relation also has the potential for devising numerical 
methods for solving BVP. Let's now consider a Brownian motion in a do- 



main D with the diffusion equation (23) and Dirichlet boundary conditions 



on dD. Now differentiate Eq. ( psj) with respect to t, and then multiply a t 
and integrate over t G (0, oo), we have 
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If one interprets P{x,t\xQ) as the probability of the Brownian particle at x 
at time t starting at xq when t = 0, then 

^dP{x, i|xo) 



u{xo) 



dx / t- 

D Jo 



dt 



-dt 



is the mean time of the particle started at xq to exit D. The left-hand-side 
of can be simplified into 

oo 

= 5{x - xo). 



dP(x,t) ^, : 
t - P{x,t) 



dt 

If now we multiply a function (f>{x), which satisfies £,(p{x) = ip{x) with (|2 
and then integrate over x £ D, we have 

dP{x,t) 



dx(l){x)S.* I t- 

D Jo 



dx{£4>{x)) / t 

D Jo 



dt 
dP{x,t 



ip{x)dx / t 

D Jo 



dt 
dP{x,t) 



dt 
-dt 



-E'- 



•xo 







dt 

tp(Xs)ds 



dt 



This shows that the rhs, in which Xt is the probabilistic Brownian motion 
and is the average along the paths of Xj started at xq, satisfies the 

inhomogeneous ODE 



£,u{x) 



-ip{x) 



This is the well-known Feynman-Kac formula 
the mean exit time of the Brownian motion X^. 



(25) 

When iIj{x) = 1, u{x) is 



4.2 Random evolution 

While a deterministic dynamic equation coupled to a white noise is called a 



stochastic differential equation and leads to Brownian motion |37], a deter- 
ministic dynamic (evolution) equation coupled to a Markov process is called 
random evolution . This is a class of stochastic models characterized by 
a system of equations like 



dP{x, 



dt 

dP{x,+) 



(F_(x)P(x, -)) - /(x)P(x, -) + gix)Pix, +) 



dt 



d_ 

dx 

— (F+(x)P(x, +)) + f{x)P{x, -) - g{x)P{x, +) 



(26) 
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There is no diffusive motion in the movement. A particle fohows determin- 
istic ODEs X = F±{x) and jumps between (+) and (— ) states. Equations 
Uke ( p6| ) have wide applications in chemistry and biology. For example, the 
stochastic averaging problem in nuclear magnetic resonance spectroscopy is 
precisely such a problem |5l|. For a recent work, see |14]. 

For large / and g, the motion is approximately x 
For extremely large / and g, the Markovian process approaches a rapidly 



gJx)F.(x)+fix)F+{x) 



varying white noise and pq ) again approaches a diffusion equation |51]. In 
spectroscopy, this corresponds to two distinct spectral lines merging into a 
single broad peak. 

For small / and g, if both -F-(x) and F^{x) have zeros, then the motion 
of the particle is still qualitatively simple: the particle will stay at a (+) fixed 
point, jump to (— ), relax to a (— ) fixed point, and stay there until jumping 
to (+) and relaxing back to the (+) fixed point, or relaxing to another (+) 
fixed point. If the two F's are arranged appropriately, the particle can be 
continuously unidirectionally transported, step by step, as demonstrated in 

One insight from this discussion is that though Eqn (^) has no appropri- 
ate stationary solution, the time-dependent solution should be well behaved. 
This points out that one should approach the time-dependent version of (pT|) 



rather than its stationary solution (Eqn. 13). 



4.3 Small diffusion and the theory of large deviations 



Is the dynamics of the degenerate equation ( pq ) the limiting behavior of 
Eqn. (^) when D± 0? This is clearly an important mathematical question 
which also has significant relevance to the modeling of muscle contractions. 
As we have stated, one way to address this question is to develop a complete 
singular perturbation theory for (|5|). There is, however, also a stochastic 
approach called theory of large deviations |13| which in recent years has 



thrown much light on the problem. A combination of both approaches 
is undoubtedly desirable. This technically very demanding task has been 



carried out in several occasions, for example in |34|. For a review, see 



5 Deterministic vs. Stochastic Motion of Molecu- 
lar Motors 

While there is no doubt that the motor protein movement is a drifted Brwo- 
nian motion, the extent of the randomness in the motion can be quantita- 
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tively characterized according to the mathematical model. Let's now again 
consider Eqn. (^) in which there is no force in the detached state (— ). As 
we have discussed above, in the limit of both D± = 0, the motion will be 
trapped at the zeros of F{x). This indicates the importance of nonzero D_ 
for the motor movement in this model, as has been repeatedly pointed out 
by Peskin et al 



On the other hand, when F-{x) 7^ 0, a motor protein can move strictly 
in one direction in a random evolution model. Diffusion plays no role in 
this mechanism. These two different modes of movement correspond nicely 
with "Brownian ratchet" and "power stroke" in the biochemical literature. 
Whether a motor protein in fact moves back-and-forth with a drift or al- 
most unidirectional consecutively can be quantitatively analyzed. Until now, 
there has no quantitative means to differentiate these two types of motion. 
The present theory offers a quantity method to address this issue. Taking 
Eqn. (H) as an example, we can introduce a function 



D 



dPjx) 
dx 



+ 



F{x)P{x) 



dx 



(27) 



as the total movement of the protein. Note that the first term is associ- 
ated with the Brownian motion and the second term is associated with the 
unidirectional movement. Hence their ratio quantitatively characterizes the 
mode of the motor movement. This integral is known as action in the theory 



of large deviations [18|. 



6 Future Work 

While the detailed mathematical analyses remain to be carried out for 
models of single motor movement, the mathematic analysis of a chain of 
motor protein is largely unknown, except for the completely rigid chain of 
motors (Huxley model). The general theory can be developed by connecting 
N motor proteins by springs. Such a "bead-and-spring" model has been the 
theoretical foundation of polymer physics The late Professor P.J. 

Flory was awarded the Nobel Prize in chemistry in 1974 for his contribution 
to this theory. The difference, however, is that a polymer is an equilibrium 
system, while a chain of motors is a "living creature". Let's denote the 
positions of motor proteins by xi, X2, xn, and the corresponding 
internal states by cJi, o"2, (Tn, where (7^ = 0,1 for the detached and 
attached state of the kth motor. We therefore have the dynamic equation 
for probability P{xi, cti, X2, (T2, xiy,aN, t)'- 
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dP 
'dt 



N 

E 

k=l I 



{(1 - a,)D^ + a,D^}^ - — + 



+ 



dxl 



+ ^ I ^(^fc-i - 2xfc + Xk+i)P 
-{(1 " <^k)f{x) +akg{x)}P{...,ak,...) 



+ {(1 - cjk)g{x) + ak)f{x)}P{..., 1 - ak, ...) 



(28) 



where ry is a spring constant. A computational analysis of such a model 
can be found in In a recent mathematical analysis, a deterministic 
counterpart of this system, a chain of bead-and-spring in a periodic force 
field, has been shown to exhibit globally phase-locked motion ||55| , 57]. Eqn. 
( p8| ) is a N-particle system which can be subjected to mean-field treatment as 



that for the N-particle Schordinger equation. As in the genesis of nonlinear 
Schrodinger equation |5^, such treatment will lead to a nonlinear diffusion 
equation ||2^, opening a possible new mathematical approach to the problem 
of muscle contraction. In connection to the theory of probability, this is 
an interacting particle system with a nonequilibrium (Gibbsian) stationary 



state [33 1, and is a natural application for the theory of large deviations |1C]. 
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